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Abstract. This lecture provides some introduction to perfect fluid dy- 
namics within the framework of general relativity. The presentation 
is based on the Carter-Lichnerowicz approach. It has the advantage 
over the more traditional approach of leading very straightforwardly to 
^ ■ important conservation laws, such as the relativistic generalizations of 

Bernoulli's theorem or Kelvin's circulation theorem. It also permits to 
, get easily first integrals of motion which are particularly useful for com- 

puting equilibrium configurations of relativistic stars in rotation or in 
binary systems. The presentation is relatively self-contained and does 
not require any a priori knowledge of general relativity. In particular, 
the three types of derivatives involved in relativistic hydrodynamics 
are introduced in detail: this concerns the Lie, exterior and covariant 
O ■ derivatives. 

cr 

1 Introduction 

> ■ 

• i-H ■ . . . . 

k>( , Relativistic fluid dynamics is an important topic of modern astrophysics in at least 
j_j ■ three contexts: (i) jets emerging at relativistic speed from the core of active galactic 
& ' nuclei or from microquasars, and certainly from gamma-ray burst central engines, 
(ii) compact stars and flows around black holes, and (iii) cosmology. Notice that 
for items (ii) and (iii) general relativity is necessary, whereas special relativity is 
sufficient for (i). 

We provide here an introduction to relativistic perfect fluid dynamics in the 
framework of general relativity, so that it is applicable to all themes (i) to (iii). 
However, we shall make a limited usage of general relativistic concepts. In par- 
ticular, we shall not use the Riemannian curvature and all the results will be 
independent of the Einstein equation. 
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We have chosen to introduce relativistic hydrodynamics via an approach devel- 
oped originally by Lichnerowicz (1941, 1955, 1967) and extended significantly by 
Carter (1973, 1979, 1989). This formulation is very elegant and permits an easy 
derivation of the relativistic generalizations of all the standard conservation laws 
of classical fluid mechanics. Despite of this, it is absent from most (all ?) text- 
books. The reason may be that the mathematical settings of Carter-Lichnerowicz 
approach is Cartan's exterior calculus, which departs from what physicists call 
"standard tensor calculus". Yet Cartan's exterior calculus is simpler than the 
"standard tensor calculus" for it does not require any specific structure on the 
spacetime manifold. In particular, it is independent of the metric tensor and its 
associated covariant derivation, not speaking about the Riemann curvature tensor. 
Moreover it is well adapted to the computation of integrals and their derivatives, 
a feature which is obviously important for hydrodynamics. 

Here we start by introducing the general relativistic spacetime as a pretty sim- 
ple mathematical structure (called manifold) on which one can define vectors and 
multilinear forms. The latter ones map vectors to real numbers, in a linear way. 
The differential forms on which Cartan's exterior calculus is based are then simply 
multilinear forms that are fully antisymmetric. We shall describe this in Sec. 2, 
where we put a special emphasis on the definition of the three kinds of derivative 
useful for hydrodynamics: the exterior derivative which acts only on differential 
forms, the Lie derivative along a given vector field and the covariant derivative 
which is associated with the metric tensor. Then in Sec. 3 we move to physics 
by introducing the notions of particle worldlinc, proper time, 4-velocity and 4- 
acceleration, as well as Lorentz factor between two observers. The hydrodynamics 
then starts in Sec. 4 where we introduce the basic object for the description of 
a fluid: a bilinear form called the stress-energy tensor. In this section, we define 
also the concept of perfect fluid and that of equation of state. The equations of 
fluid motion are then deduced from the local conservation of energy and momen- 
tum in Sec. 5. They are given there in the standard form which is essentially 
a relativistic version of Euler equation. From this standard form, we derive the 
Carter-Lichnerowicz equation of motion in Sec. 6, before specializing it to the case 
of an equation of state which depends on two parameters: the baryon number 
density and the entropy density. We also show that the Newtonian limit of the 
Carter-Lichnerowicz equation is a well known alternative form of the Euler equa- 
tion, namely the Crocco equation. The power of the Carter-Lichnerowicz approach 
appears in Sec. 7 where we realize how easy it is to derive conservation laws from it, 
among which the relativistic version of the classical Bernoulli theorem and Kelvin's 
circulation theorem. We also show that some of these conservation laws are useful 
for getting numerical solutions for rotating relativistic stars or relativistic binary 
systems. 

2 Fields and derivatives in spacetime 

It is not the aim of this lecture to provide an introduction to general relativity. 
For this purpose we refer the reader to two excellent introductory textbooks which 
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have recently appeared: (Hartle 2003) and (Carrol 2004). Here we recall only 
some basic geometrical concepts which are fundamental to a good understanding 
of relativistic hydrodynamics. In particular we focus on the various notions of 
derivative on spacetime. 

2.1 The spacetime of general relativity 

Relativity has performed the fusion of space and time, two notions which were 
completely distinct in Newtonian mechanics. This gave rise to the concept of 
spacetime, on which both the special and general theory of relativity are based. 
Although this is not particularly fruitful (except for contrasting with the rela- 
tivistic case), one may also speak of spacetime in the Newtonian framework. The 
Newtonian spacetime Ai is then nothing but the affine space R 4 , foliated by the 
hyperplanes E t of constant absolute time t: these hyperplanes represent the ordi- 
nary 3-dimensional space at successive instants. The foliation (£ t ) tG R is a basic 
structure of the Newtonian spacetime and does not depend upon any observer. 
The worldline £ of a particle is the curve in M generated by the successive posi- 
tions of the particle. At any point A E C, the time read on a clock moving along 
C is simply the parameter t of the hyperplanc S t that intersects C at A. 

The spacetime M of special relativity is the same mathematical space as the 
Newtonian one, i.e. the affine space R 4 . The major difference with the Newtonian 
case is that there does not exist any privileged foliation (X t ) tG R. Physically this 
means that the notion of absolute time is absent in special relativity. However 
M is still endowed with some absolute structure: the metric tensor g and the 
associated light cones. The metric tensor is a symmetric bilinear form g on Ai, 
which defines the scalar product of vectors. The null (isotropic) directions of g 
give the worldlines of photons (the light cones) . Therefore these worldlincs depend 
only on the absolute structure g and not, for instance, on the observer who emits 
the photon. 

The spacetime M of general relativity differs from both Newtonian and special 
relativistic spacetimes, in so far as it is no longer the affine space R 4 but a more 
general mathematical structure, namely a manifold. A manifold of dimension 4 
is a topological space such that around each point there exists a neighbourhood 
which is homeomorphic to an open subset of R 4 . This simply means that, locally, 
one can label the points of M in a continuous way by 4 real numbers (a; Q ) a e{o.i,2,3} 
(which are called coordinates). To cover the full M, several different coordinates 
patches (charts in mathematical jargon) can be required. 

Within the manifold structure the definition of vectors is not as trivial as within 
the affine structure of the Newtonian and special relativistic spacetimes. Indeed, 
only infinitesimal vectors connecting two infinitely close points can be defined a 
priori on a manifold. At a given point p E M, the set of such vectors generates a 
4-dimensional vector space, which is called the tangent space at the point p and 
is denoted by T p (Ai). The situation is therefore different from the Newtonian or 
special relativistic one, for which the very definition of an affine space provides 
a unique global vector space. On a manifold there are as many vector spaces as 
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Fig. 1. The vectors at two points p and q on the spacetime manifold M belong to two 
different vector spaces: the tangent spaces T P (M) and T q (M). 



points p (cf. Fig. 1). 

Given a vector basis (e a ) ae r 0i i,2,3} of T p (M), the components of a vector 
v 6 T P (M) on this basis are denoted by v a : v — v a e a , where we have employed 
Einstein's convention for summation on repeated indices. It happens frequently 
that the vector basis is associated to a coordinate system (x a ) on the manifold, in 
the following way. If dx a is the (infinitesimal) difference of coordinates between p 
and a neighbouring point q, the components of the vector pq' with respect to the 
basis (e a ) are exactly dx a . The basis that fulfills this property is unique and is 
called the natural basis associated with the coordinate system (x a ). It is usually 
denoted by (d/dx a ), which is a reminiscence of the intrinsic definition of vectors 
on a manifold as differential operators acting on scalar fields. 

As for special relativity, the absolute structure given on the spacetime manifold 
M. of general relativity is the metric tensor g. It is now a field on M.: at each 
point peM, g{p) is a symmetric bilinear form acting on vectors in the tangent 
space T P (M): 



It is demanded that the bilinear form g is not degenerate and is of signature 
(— , +, +, +). It thus defines a scalar product on T p (A4), which justifies the notation 
u ■ v for g(u,v). The isotropic directions of g give the local light cones: a vector 
v G T p {M) is tangent to a light cone and called a null or lightlike vector iff v-v = 0. 
Otherwise, the vector is said to be timelike iff v ■ v < and spacelike iff v ■ v > 0. 



g(p): T p (M)xT p (M) 
(u, v) 



R 

g(u,v) 



=: u ■ v. 



(2.1) 
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2.2 Tensors 

Let us recall that a linear form at a given point p <E M is an application 

" T >M - ? A „ (2.2) 
t> i — > {u>, v) := U3{v), v ' 

that is linear. The set of all linear forms at p forms a vector space of dimension 4, 
which is denoted by T p (M)* and is called the dual of the tangent space T p {M). In 
relativistic physics, an abundant use is made of linear forms and their generaliza- 
tions: the tensors. A tensor of type (k,£), also called tensor k times contravariant 
and £ times covariant, is an application 

T : T P {M)* x • • • x T P {M)* x T P (M) x • • • x T P (M) — > R 

k times £ times 

(u>i,...,u> fe , vi,..., v e ) i — > r(wi,...,w fe , 

Vl, •••,«<) 

(2.3) 

that is linear with respect to each of its arguments. The integer k + I is called 
the valence of the tensor. Let us recall the canonical duality T p (A4)** = T p (A4), 
which means that every vector v can be considered as a linear form on the space 
T p (Ad)*, defining the application v : T p (M)* — > R, lj (w, v). Accordingly a 
vector is a tensor of type (1,0). A linear form is a tensor of type (0, 1) and the 
metric tensor g is a tensor of type (0, 2). 

Let us consider a vector basis of T p {M) 1 (e Q ), which can be either a natural 
basis (i.e. related to some coordinate system) or not (this is often the case for bases 
orthonormal with respect to the metric g). There exists then a unique quadruplet 
of 1-forms, (e Q ), that constitutes a basis of the dual space T p {M)*~ and that verifies 

(e a ,e fi )=6 a p , (2.4) 

where 5 a l3 is the Kronecker symbol. Then we can expand any tensor T of type 
(k,£) as 

! " " ' (2-5) 



T = T ai - ak pi fjf e ai ® . . . <g> e ak ® ® . . . eg) 



where the tensor product e ai ® . . . ® e Qfc ® e® 1 ® . . . ® e^ e is the tensor of type 
(k, £) for which the image of (u>i, . . . ,u>k,vi, ■ ■ ■ ,ve) as in (2.3) is the real number 



k 



n<^,e ai ) x JJ(e*,tf,-). (2.6) 

»=i j=i 

Notice that all the products in the above formula are simply products in R. The 
4 fe+£ scalar coefficients j ,ai --- a: '= i n (2.5) are called the components of the tensor 
T with respect to the basis (e a ), or with respect to the coordinates (x a ) if (e Q ) 
is the natural basis associated with these coordinates. These components are 
unique and fully characterize the tensor T. Actually, in many studies, a basis is 
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assumed (mostly a natural basis) and the tensors are always represented by their 
components. This way of presenting things is called the index notation, or the 
abstract index notation if the basis is not specified (e.g. Wald 1984). We shall not 
use it here, sticking to what is called the index-free notation and which is much 
better adapted to exterior calculus and Lie derivatives. 

The notation v a already introduced for the components of a vector v is of 
course the particular case (k — 1,£ = 0) of the general definition given above. 
For a linear form u>, the components uo a are such that u> = uo a e a [Eq. (2.5) with 
(k = 0,£=l)]. Then 

(u,,v)=u> a v a . (2.7) 

Similarly the components g a $ of the metric tensor g are defined by g — g a p e"(g>e^ 
[Eq. (2.5) with (k = 0, £ = 2)] and the scalar products are expressed in terms of 
the components as 

g(u,v) =g al3 u a vf 3 . (2.8) 

2.3 Scalar fields and their gradients 

A scalar field on the spacetime manifold M is an application / : M — > M. If / 
is smooth, it gives rise to a field of linear forms (such fields are called 1-forms), 
called the gradient of f and denoted df. It is defined so that the variation of / 
between two neighbouring points p and q is 1 

<V = f(q)-m = (df,m- (2-9) 

Let us note that, in non-relativistic physics, the gradient is very often considered as 
a vector and not as a 1-form. This is because one associates implicitly a vector Q to 
any 1-form oj thanks to the Euclidean scalar product of M 3 , via V)/eR 3 , (oj, v) = 
Q ■ v. Accordingly, the formula (2.9) is rewritten as df — V/ • p§. But one shall 
keep in mind that, fundamentally, the gradient is a 1-form and not a vector. 

If (x a ) is a coordinate system on M and (e a — d/dx a ) the associated natural 
basis, then the dual basis is constituted by the gradients of the four coordinates: 
e a = dx a . The components of the gradient of any scalar field / in this basis are 
then nothing but the partial derivatives of /: 

d/ = (d/) Q d.T Q with (df) a = ^L. (2.10) 

2.4 Comparing vectors and tensors at different spacetime points: various 
derivatives on hA 

A basic concept for hydrodynamics is of course that of vector field. On the manifold 
A4, this means the choice of a vector v(p) in T p (M) for each p 6 M. We denote by 



1 do not confuse the increment df of / with the gradient 1-form df: the boldface d is used to 
distinguish the latter from the former 



Eric Gourgoulhon: An introduction to relativistic hydrodynamics 7 




Fig. 2. Geometrical construction of the Lie derivative of a vector field: given a small 
parameter A, each extremity of the arrow Xv is dragged by some small parameter e along 
u, to form the vector denoted by <!> e (A?;). The latter is then compared with the actual 
value of Xv at the point q, the difference (divided by Ae) defining the Lie derivative Ca v. 

T{M) the space of all smooth vector fields on Ai 2 . The derivative of the vector 
field is to be constructed for the variation Sv of v between two neighbouring 
points p and q. Naively, one would write Sv = v(q) — v(p), as in (2.9). However 
v(q) and v(p) belong to different vector spaces: T q {M) and T p {M) (cf. Fig. 1). 
Consequently the subtraction v(q) — v(p) is ill defined, contrary of the subtraction 
of two real numbers in (2.9). To proceed in the definition of the derivative of a 
vector field, one must introduce some extra-structure on the manifold A4: this can 
be either another vector field u, leading to the derivative of v along u which is 
called the Lie derivative, or a connection V (usually related to the metric tensor 
g) , leading to the covariant derivative ViJ. These two types of derivative generalize 
straightforwardly to any kind of tensor field. For the specific kind of tensor fields 
constituted by differential forms, there exists a third type of derivative, which does 
not require any extra structure on A4 : the exterior derivative. We will discuss the 
latter in Sec. 2.5. In the current section, we shall review successively the Lie and 
covariant derivatives. 

2.4.1 Lie derivative 

The Lie derivative is a very natural operation in the context of fluid mechanics. 
Indeed, consider a vector field u on Ai, called hereafter the flow. Let v be another 
vector field on Ai, the variation of which is to be studied. We can use the flow u 
to transport the vector v from one point p to a neighbouring one q and then define 



2 The experienced reader is warned that T(M) does not stand for the tangent bundle of M 
(it rather corresponds to the space of smooth cross-sections of that bundle). No confusion may 
arise since we shall not use the notion of bundle. 
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rigorously the variation of v as the difference between the actual value of v at q 
and the transported value via u. More precisely the definition of the Lie derivative 
of v with respect to u is as follows (see Fig. 2). We first define the image $ £ (p) of 
the point p by the transport by an infinitesimal "distance" e along the field lines 
of u as $ e (p) = q, where q is the point close to p such that p$ = eu(p). Besides, 
if we multiply the vector v(p) by some infinitesimal parameter A, it becomes an 
infinitesimal vector at p. Then there exists a unique point p' close to p such that 
Xv(p) = pp. We may transport the point p' to a point q' along the field lines of u 
by the same "distance" e as that used to transport p to q: q' = <& £ (p') (see Fig. 2). 
qq' is then an infinitesimal vector at q and we define the transport by the distance 
e of the vector v(p) along the field lines of u according to 

$ e (v(p)):=jM'. (2.11) 

$ e (v(p)) is vector in T q {M). We may then subtract it from the actual value of 
the field v at q and define the Lie derivative of v along u by 

C u v := lim - [v(q) - $ e («(p))] . (2.12) 

If we consider a coordinate system (x a ) adapted to the field u in the sense 
that u = eb where e is the first vector of the natural basis associated with the 
coordinates (x a ), then the Lie derivative is simply given by the partial derivative 
of the vector components with respect to x°: 

tW = &S- (2 - 13) 
In an arbitrary coordinate system, this formula is generalized to 

£ a v"=u^-v»^-, (2.14) 
dx» dx» ' 

where use has been made of the standard notation Cn v a := (Ca v) a . 

The Lie derivative is extended to any tensor field by (i) demanding that for a 
scalar field /, Ca f = (df,u) and (ii) using the Leibniz rule. As a result, the Lie 
derivative Ca T of a tensor field T of type (k, £) is a tensor field of the same type, 
the components of which with respect to a given coordinate system (x a ) are 

h i 



f rpai...a k _ ai u r r a 1 ...a k _ rpai--- a...a k 

L - ul Pi...0t ~ Pi-fa Z^ 1 01 



i—1 

f 

+E rai - a v... T ... A tosr- ( 2 - 15 ) 



In particular, for a 1-form, 



(2.16) 
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2.4.2 Covariant derivative 

The variation Sv of the vector field v between two neighbouring points p and q 
can be defined if some affine connection is given on the manifold M. The latter 
is an operator 

V : T(M) x T(M) — > T(M) 
(u,v) i — > V a v 

that satisfies all the properties of a derivative operator (Leibniz rule, etc.), which 
we shall not list here (see e.g. Wald 1984). The variation of v (with respect to the 
connection V) between two neighbouring points p and q is then defined by 

5v:=V n v. (2.18) 

One says that v is transported parallelly to itself between p and q iff 5v = 0. 
From the manifold structure alone, there exists an infinite number of possible 
connections and none is preferred. Taking account the metric tensor g changes 
the situation: there exists a unique connection, called the Levi-Civita connection, 
such that the tangent vectors to the geodesies with respect to g are transported 
parallelly to themselves along the geodesies. In what follows, we will make use 
only of the Levi-Civita connection. 

Given a vector field v and a point p e M, we can consider the type (1,1) 
tensor at p denoted by V'u(p) and defined by 

W(p) : T P (M)* xT p (M) — R . . 

(«,«) ^ («, (V ac v)(p)) ' (Z - W) 

where u c is a vector field that performs some extension of the vector u in the 
neighbourhood of p: u c (p) = u. It can be shown that the map (2.19) is inde- 
pendent of the choice of u c . Therefore V«(p) is a type (1, 1) tensor at p which 
depends only on the vector field v. By varying p we get a type (1,1) tensor field 
denoted and called the covariant derivative of v. 

As for the Lie derivative, the covariant derivative is extended to any tensor 
field by (i) demanding that for a scalar field V/ = d/ and (ii) using the Leibniz 
rule. As a result, the covariant derivative of a tensor field T of type (k,£) is a 
tensor field VT of type (k,£+ 1). Its components with respect a given coordinate 
system (x a ) are denoted 

V 7 T--<V.A : = (W) ai - a V,.. ft7 (2.20) 
(notice the position of the index 7 !) and are given by 

k i 



vy rpcei...ctk _ u r r a 1 ...a k , \ ^ poj rp' 

y f 1 - Q xl L 01-01 + 2-^ 1 I" 1 01-01 



7ft J 0i...cr...0 e i 
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where the coefficients T" ^ are the Christoffel symbols of the metric g with respect 
to the coordinates {x a ). They are expressible in terms of the partial derivatives 
of the components of the metric tensor, via 

A distinctive feature of the Levi-Civita connection is that 



Vg = . (2.23) 



Given a vector field u and a tensor field T of type (k,£), we define the covariant 
derivative ofT along it as the generalization of (2.17): 

V fi T:=VT(_,S). (2.24) 

k-\-£ slots 

Notice that T is a tensor of the same type (k, £) as T and that its components 
are 

(VuTr-<V.. A = ^n-^..^. (2.25) 
2.5 Differential forms and exterior derivatives 

The differential forms or n-forms are type (0, n) tensor fields that are antisym- 
metric in all their arguments. Otherwise stating, at each point p e they 
constitute antisymmetric multilinear forms on the vector space T p (A4)). They 
play a special role in the theory of integration on a manifold. Indeed the primary 
definition of an integral over a manifold of dimension n is the integral of a n-form. 
The 4-dimensional volume element associated with the metric g is a 4-form, called 
the Levi-Civita alternating tensor. Regarding physics, it is well known that the 
electromagnetic field is fundamentally a 2-form (the Faraday tensor F); besides, 
we shall see later that the vorticity of a fluid is described by a 2-form, which plays 
a key role in the Carter-Lichnerowicz formulation. 

Being tensor fields, the n-forms are subject to the Lie and covariant derivations 
discussed above. But, in addition, they are subject to a third type of derivation, 
called exterior derivation. The exterior derivative of a n-form u; is a (n + l)-form 
which is denoted do;. In terms of components with respect to a given coordinate 
system (x a ), da; is defined by 

0-form (scalar field) : (dw) Q = (2.26) 

l-fo„n : «M„, = fe-|j (2.2T) 
2 , orm : (d ^, = ^ + ^ + ^ „ 28) 
etc... (2.29) 
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It can be easily checked that these formulas, although expressed in terms of partial 
derivatives of components in a coordinate system, do define tensor fields. More- 
over, the result is clearly antisymmetric (assuming that u) is), so that we end up 
with (n + l)-forms. Notice that for a scalar field (0-form), the exterior derivative 
is nothing but the gradient 1-form df already defined in Sec. 2.3. Notice also that 
the definition of the exterior derivative appeals only to the manifold structure. It 
does not depend upon the metric tensor g, nor upon any other extra structure on 
M. We may also notice that all partial derivatives in the formulas (2.26)-(2.28) 
can be replaced by covariant derivatives (thanks to the symmetry of the Christoffcl 
symbols). 

A fundamental property of the exterior derivation is to be nilpotent: 



dda; = . (2.30) 



A n-form u> is said to be closed iff du> — 0, and exact iff there exists a (n — 1)- 
form er such that u> = der. From property (2.30), an exact n-form is closed. The 
Poincare lemma states that the converse is true, at least locally. 

The exterior derivative enters in the well known Stokes' theorem: if I? is a 
submanifold of M of dimension d (d e {1, 2,3,4}) that has a boundary (denoted 
dV), then for any (d — l)-form u>, 



u>= du>. (2.31) 

dV JV 

Note that dV is a manifold of dimension d — 1 and du; is a rf-form, so that each 
side of (2.31) is (of course !) a well defined quantity, as the integral of a n-form 
over a n-dimcnsional manifold. 

Another very important formula where the exterior derivative enters is the 
Cartan identity, which states that the Lie derivative of a n-form uj along a vector 
field u is expressible as 



= u ■ da; + d(u ■ u>) . (2.32) 



In the above formula, a dot denotes the contraction on adjacent indices, i.e. u ■ u) 
is the (n — l)-form u)(v, .,...,.), with the n — 1 last slots remaining free. Notice 
that in the case where a; is a 1-form, Eq. (2.32) is readily obtained by combining 
Eqs. (2.16) and (2.27). In this lecture, we shall make an extensive use of the 
Cartan identity. 



3 Worldlines in spacetime 

3.1 Proper time, 4-velocity and 4-acceleration 

A particle or "point mass" is fully described by its mass m > and its worldlinc 
£ in spacetime. The latter is postulated to be timelike, i.e. such that any tangent 
vector is timelike. This means that C lies always inside the light cone (see Fig. 3). 
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Fig. 3. Worldline C of a particle, with the 4-velocity u = c 1 dx/dr and the 4-acceleration 
a. 



The proper time dr corresponding to an elementary displacement 3 dx along C is 
nothing but the length, as given by the metric tensor, of the vector dx (up to a c 
factor) : 

cdr = \/—g(dx, dx). (3.1) 
The ^-velocity of the particle is then the vector defined by 



1 dx 

c dr 



(3.2) 



By construction, u is a vector tangent to the worldline C and is a unit vector with 
respect to the metric g: 

u ■ u = — 1. (3-3) 

Actually, u can be characterized as the unique unit tangent vector to C oriented 
toward the future. Let us stress that the 4-velocity is intrinsic to the particle under 
consideration: contrary to the "ordinary" velocity, it is not defined relatively to 
some observer. 

The 4~ acceleration of the particle is the covariant derivative of the 4-velocity 
along itself: 

' ' (3-4) 



V,7 W 



Since it is a unit vector, it follows that 

u ■ a = 0, 



(3.5) 



i.e. a, is orthogonal to u with respect to the metric g (cf. Fig. 3). In particular, a 
is a spacelike vector. Again, the 4-acceleration is not relative to any observer, but 
is intrinsic to the particle. 



3 we denote by dx the infinitesimal vector between two neighbouring points on C, but it should 
be clear that this vector is independent of any coordinate system (x a ) on M. 
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Fig. 4. Relative velocity U = dx/dro of a particle of 4- velocity u with respect to an 
observer of 4- velocity uo- U enters in the orthogonal decomposition of u with respect to 
uq, via u — r(ito + c _1 U). NB: contrary to what the figure might suggest, dro > dr. 



3.2 Observers, Lorentz factors and relative velocities 

Let us consider an observer Oq (treated as a point mass particle) of worldline Cq. 
Let us recall that, following Einstein's convention for the definition of simultaneity, 
the set of events that are considered by Oq as being simultaneous to a given event 
p on his worldline is a hypersurface of M. which is orthogonal (with respect to g) 
to Cq at p. 

Let O be another observer, whose worldline C intersects that of Oq at p. Let us 
denote by tq (resp. r) the proper time of Co (resp. O). After some infinitesimal 
proper time dr, O is located at the point q (cf. Fig. 4). Let then to + dro be 
the date attributed by Co to the event q according to the simultaneity convention 
recalled above. The relation between the proper time intervals dro and dr is 

dro = Tdr, (3-6) 

where T is the Lorentz factor between the observers Co and O. We can express T 
in terms of the 4- velocities Uq and u of Oo and O. Indeed, let dx the infinitesimal 
vector that is orthogonal to «o and links Cq to q (cf. Fig. 4). Since ilo and u are 
unit vectors, the following equality holds: 

cdru = cdTQUQ + dx. (3-7) 

Taking the scalar product with Uq, and using (3.6) as well as Uq ■ dx = results 

= —mq • u . (3.8) 



Hence from a geometrical point of view, the Lorentz factor is nothing but (minus) 
the scalar product of the unit vectors tangent to the two observers' worldlincs. 



14 



Title : will be set by the publisher 



The velocity of O relative to Oo is simply the displacement vector dx divided 
by the elapsed proper time of O , dr : 

dx 



U := 



dr ' 



(3.9) 



U is the "ordinary" velocity, by opposition to the 4-velocity u. Contrary to the 
latter, which is intrinsic to O, U depends upon the observer Oq. Geometrically, 
U can be viewed as the part of u that is orthogonal to u n , since by combining 
(3.6) and (3.7), we get 



a = r \u + -u 



with u ■ U = 0. 



(3.10) 



Notice that Eq. (3.8) is a mere consequence of the above relation. The scalar square 
of Eq. (3.10), along with the normalization relations u u = —1 and u • «o = —1, 
leads to 



r = ( i - \u ■ u 



-1/2 



which is identical to the well-known expression from special relativity. 



(3.11) 



4 Fluid stress-energy tensor 



4.1 General definition of the stress-energy tensor 

The stress-energy tensor T is a tensor field on M which describes the matter 
content of spacetime, or more precisely the energy and momentum of matter, at 
a macroscopic level. T is a tensor field of type (0, 2) that is symmetric (this 
means that at each point p € Ai, T is a symmetric bilinear form on the vector 
space T p {M)) and that fulfills the following properties: given an observer O of 
4-velocity Uo, 

• the matter energy density as measured by Oq is 

E = T(u ,u ); (4.1) 

• the matter momentum density as measured by Co is 

P= --T(u ,ei)ei, (4.2) 

c 

where (Si) is an orthonormal basis of the hyperplanc orthogonal to Oo's 
worldline (rest frame of Co); 

• the matter stress tensor as measured by Oo is 

Sij =T(e i ,e j ), (4.3) 

i.e. T(e>i, ej) is the force in the direction acting on the unit surface whose 
normal is e, . 
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4.2 Perfect fluid stress-energy tensor 

The perfect fluid model of matter relies on a field of 4- velocities u, giving at each 
point the 4-velocity of a fluid particle. Moreover the perfect fluid is characterized 
by an isotropic pressure in the fluid frame (i.e. Sij = pSij for the observer whose 
4-velocity is u). More precisely, the perfect fluid model is entirely defined by the 
following stress-energy tensor: 



T = (pc 2 +p)u®u+pg , (4.4) 



where p and p are two scalar fields, representing respectively the matter energy 
density (divided by c 2 ) and the pressure, both measured in the fluid frame, and u 
is the 1-form associated to the 4-velocity u by the metric tensor g: 

* 1 TpiM) ^ V ^ - - (4-5) 
v i — ► g(u,v) = u-v. 

In terms of components with respect to a given basis (S a ), if u — u a S a and 
if (e a ) is the 1-form basis dual to (S a ) (cf. Sec. 2.2), then u = u a e a , with 
u a = 9ap uli ■ In Eq. (4.4), the tensor product u®u stands for the bilinear form 
(v,w) i ► (u, v) (u, w) = (u ■ v)(u ■ w) [cf. (2.6)] . 

According to Eq. (4.1) the fluid energy density as measured by an observer 
Oq of 4-velocity Uq is E = T(u n ,u ) = (pc 2 + p)(u ■ u ) 2 + pg(u ,uo)- Since 
u ■ Hq = — r, where T is the Lorentz factor between the fluid and Oq [Eq. (3.8)], 
and gr(tto,tto) = —1, we get 

E = T 2 (pc 2 +p)-p. (4.6) 

The reader familiar with the formula E = Tmc 2 may be puzzled by the T 2 factor 
in (4.6). However he should remind that E is not an energy, but an energy per 
unit volume: the extra T factor arises from "length contraction" in the direction 
of motion. 

Similarly, by applying formula (4.2), we get the fluid momentum density as 
measured by the observer O : cp = —T(u , e^) = — [(pc 2 + p)(u ■ Uo)(u ■ Si) + 
pg(uo, e"i)] ei, with u ■ u = —T, g(uo, e-i) = and (u ■ Si) Si being the projection 
of u orthogonal to tto: according to (3.10), (u ■ Si) Si = T/cU, where U is the 
fluid velocity relative to Co- Hence 

P = r2 {p + ^)u- (4-7) 

Finally, by applying formula (4.3), we get the stress tensor as measured by 
the observer Oq: Sij = T(Si,Sj) = (pc 2 + p)(u ■ Si)(u ■ Sj) + pg(Si,Sj), with 
u-Si — T l 'cSi-U [thanks to Eq. (3.10) and Si ■ tt = 0] andg^e^e,) =Sij. Hence 

S^^pS^+T 2 ( p+ ^WU\ (4.8) 
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where U l is the i-th component of the velocity U with respect to the orthonormal 
triad (e^): U = U 1 ^. 

Notice that if the observer O is comoving with the fluid, then u = u, T = 1, 
U = and Eqs. (4.6), (4.7) and (4.8) reduce to 

E = pc 2 , p = 0, S tJ =pS l3 . (4.9) 

We thus recover the interpretation of the scalar fields p and p given above. 



4.3 Concept of equation of state 

Let us assume that at the microscopic level, the perfect fluid is constituted by N 
species of particles (N > 1), so that the energy density pc 2 is a function of the 
number densities ua of particles of species A (A e {1,2,..., TV}) in the fluid rest 
frame (proper number density) and of the entropy density s in the fluid rest frame: 



pc 2 = e(s,ni,n 2 , . . . ,n N ) . (4-10) 



The function e is called the equation of state (EOS) of the fluid. Notice that e 
is the total energy density, including the rest-mass energy: denoting by m A the 
individual mass of particles of species A, we may write 

£ = ^m A n A c 2 +£i„t, (4.11) 

A 

where £i nt is the "internal" energy density, containing the microscopic kinetic 
energy of the particles and the potential energy resulting from the interactions 
between the particles. 

The first law of thermodynamics in a fixed small comobile volume V writes 

d£ = TdS + ^2n A dN A , (4.12) 

A 

where 

• £ is the total energy in volume V: £ = eV, 

• S is the total entropy in V: S = sV, 

• Na is the number of particles of species A in V: Na = uaV, 

• T is the thermodynamical temperature, 

• p A is the relativistic chemical potential of particles of particles of species A; 
it differs from the standard (non-relativistic) chemical potential fi A by the 
mass m A : p A = p A + m A , reflecting the fact that £ includes the rest-mass 
energy. 
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Replacing £ , S and Na by their expression in terms of e, s, and V leads 
to d{eV) = Td(sV) + *£ A (j, A d{n A V). Since V is held fixed, the first law of 
thermodynamics becomes 



de = Tds + ^2fi A dn A 



(4.13) 



Consequently, T and \i A can be expressed as partial derivatives of the equation of 
state (4.10): 



.' de 
~ds 



and 



(4.14) 



Actually these relations can be taken as definitions for the temperature T and 
chemical potential [i A . This then leads to the relation (4.13), which we will call 
hereafter the first law of thermodynamics. 



5 Conservation of energy and momentum 

5.1 General form 

We shall take for the basis of our presentation of relativistic hydrodynamics the 
law of local conservation of energy and momentum of the fluid, which is assumed 
to be isolated: 



V ■ T = | . (5.1) 



V • T stands for the covariant divergence of the fluid stress-energy tensor T. This 
is a 1-form, the components of which in a given basis are 

( V • 1% = V^ Q = g» v V v T m . (5.2) 

For a self-gravitating fluid, Eq. (5.1) is actually a consequence of the fundamental 
equation of general relativity, namely the Einstein equation. Indeed the latter re- 
lates the curvature associated with the metric g to the matter content of spacetime, 
according to 

G~T, (5.3) 

where G is the so-called Einstein tensor, which represents some part of the Ric- 
mann curvature tensor of (M,g). A basic property of the Einstein tensor is 

V • G = (this follows from the so-called Bianchi identities, which are pure 
geometric identities regarding the Riemann curvature tensor). Thus it is imme- 
diate that the Einstein equation (5.3) implies the energy-momentum conservation 
equation (5.1). Note that in this respect the situation is different from that of 
Newtonian theory, for which the gravitation law (Poisson equation) does not im- 
ply the conservation of energy and momentum of the matter source. We refer the 
reader to § 22.2 of (Hartle 2003) for a more extended discussion of Eq. (5.1), in 
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particular of the fact that it corresponds only to a local conservation of energy 
and momentum. 

Let us mention that there exist formulations of relativistic hydrodynamics that 
do not use (5.1) as a starting point, but rather a variational principle. These 
Hamiltonian formulations have been pioneered by Taub (1954) and developed, 
among others, by Carter (1973, 1979, 1989), as well as Comer & Langlois (1993). 

5.2 Application to a perfect fluid 

Substituting the perfect fluid stress-energy tensor (4.4) in the energy-momentum 
conservation equation (5.1), and making use of (2.23) results in 

V-T = ^> [V c (e +p) + (e +p)V • u]u+ (e + p)a + Vp = 0, (5.4) 

where a is the 1-form associated by the metric duality [cf. Eq. (4.5) with u 
replaced by a] to the fluid 4-acceleration a, = u [Eq. (3.4)]. The scalar V • u 
is the covariant divergence of the 4- velocity vector: it is the trace of the covariant 
derivative Vu, the latter being a type (1, 1) tensor: V • u = V cr u a . Notice that 
Vp in Eq. (5.4) is nothing but the gradient of the pressure field: Vp = dp (cf. 
item (i) in Sec. 2.4.2). 

5.3 Projection along u 

Equation (5.4) is an identity involving 1-forms. If we apply it to the vector u, we 
get a scalar field. Taking into account {u,u) = u u = — 1 and (o, u) = a ■ u = 
[Eq. (3.5)], the scalar equation becomes 

(V-T,«}=0 V a e = -(e+p)V ■ u. (5.5) 

Notice that e = (de,u) — C^e (cf. item (i) at the end of Sec. 2.4.1). Now, 
the first law of thermodynamics (4.13) yields 

V u e = TVns + ^2^ A Vnn Al (5.6) 

A 

so that Eq. (5.5) can be written as 

(V-T,m) = TV • (su) + M A V • (n A u) 

A 

+ ^e+p-Ts-^2(i A n A ^jv ■ u = 0. (5.7) 

Now, we recognize in Q := e+p—Ts the free enthalpy (also called Gibbs free energy) 
per unit volume. It is well known that the free enthalpy G = QV = E + PV — TS 
(where V is some small volume element) obeys the thermodynamic identity 

G = J2^ A N A , (5.8) 

A 
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from which we get Q = ^2 a ^ An A, i.e. 



p = Ts + h A ua 



(5.9) 



This relation shows that p is a function of (s, n 1; . . . , n^) which is fully determined 
by e(s, ni, . . . , njv) [recall that T and [i A are nothing but partial derivatives of the 
latter, Eq. (4.14)]. Another way to get the identity (5.9) is to start from the first 
law of thermodynamics in the form (4.12), but allowing for the volume V to vary, 
i.e. adding the term — pdV to it: 



d£ = T dS - p dV + V A dN A . 



(5.10) 



Substituting £ — eV, S — sV and Na = tiaV is this formula and using (4.13) 
leads to (5.9). 

For our purpose the major consequence of the thermodynamic identity (5.9) is 
that Eq. (5.7) simplifies substantially: 



(V-T,t*> 







TV • (su) + V A V ' ( n Au) = 



(5.11) 



In this equation, cV • (su) is the entropy creation rate (entropy created per unit 
volume and unit time in the fluid frame) and cV • (uau) is the particle creation 
rate of species A (number of particles created par unit volume and unit time in 
the fluid frame). This follows from 



cV • (n A u) = c(VariA + n A V ■ u) 



dn A , 1 dV 1 d(n A V) 

-^ + nA v^ = v^^' (5 - 12) 



where r is the fluid proper time and where we have used the expansion rate formula 

(5.13) 



„ - 1 dV 



V being a small volume element dragged along by u. Equation (5.11) means that 
in a perfect fluid, the only process that may increase the entropy is the creation 
of particles. 



5.4 Projection orthogonally to u: relativistic Euler equation 

Let us now consider the projection of (5.4) orthogonally to the 4-velocity. The 
projector orthogonal to u is the operator P := 1 + u <g> u: 



T P (M) 



T P (M) 

v + (u ■ v)u. 



(5.14) 
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Combining P to the 1-form (5.4), and using uoP 
to the 1-form equation 



as well as a o P = a, leads 



(V • T) o P = 



(e + p)a = -Vp- (Vu p)u 



(5.15) 



This is clearly an equation of the type " ma = F ", although the gravitational 
"force" is hidden in the covariant derivative in the derivation of a from u. We 
may consider that (5.15) is a relativistic version of the classical Euler equation. 

Most textbooks stop at this point, considering that (5.15) is a nice equation. 
However, as stated in the Introduction, there exists an alternative form for the 
equation of motion of a perfect fluid, which turns out to be much more useful 
than (5.15), especially regarding the derivation of conservation laws: it is the 
Carter-Lichnerowicz form, to which the rest of this lecture is devoted. 



6 Carter-Lichnerowicz equation of motion 



6.1 Derivation 

In the right hand-side of the relativistic Euler equation (5.15) appears the gradient 
of the pressure field: Vp = dp. Now, by deriving the thermodynamic identity (5.9) 
and combining with the first law (4.13), we get the relation 



(6.1) 



which is known as the Gibbs-Duhem relation. We may use this relation to express 
Vp in terms of VT and V[i A in Eq. (5.15). Also, by making use of (5.9), we may 
replace e + p by Ts + ^2 A pL A nA- Hence Eq. (5.15) becomes 

[ts + H^aJ a = -sVT-J2 n A V^ A - (^V a T + ^ n A V a ^ u. (6.2) 



dp = sdT + ^ riAd^ 



Writing a = V^tt and reorganizing slightly yields 

s [Va(Tu) + VT] + J2 n A [V n (p A u) + V M A ] = 0. (6.3) 

A 

The next step amounts to noticing that 

Va(Tu) =Ca(Tu). (6.4) 

This is easy to establish, starting from expression (2.16) for the Lie derivative of 
a 1-form, in which we may replace the partial derivatives by covariant derivatives 
[thanks to the symmetry of the Christoffel symbols, cf. Eq. (2.21)]: 

jC-u (Tu a ) = u^V„(Tu a ) + Tu^ aU ^ (6.5) 
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i.e. 



Al (Tu) = V u (Tu) +Tu- Vu. 



Now, from u ■ u = —1, we get u ■ Vtt = 0, which establishes (6.4). 
On the other side, the Cartan identity (2.32) yields 

£>a (Tu) = u ■ d(Tu) + d[T (u, u)]=u- d(Tu) - dT. 

Combining this relation with (6.4) (noticing that dT = VT), we get 

V a (Tu) + VT = u-d(Tu). 

Similarly, 

Va(p A u) + V[i A = u ■ d(/j, A u). 



(6.6) 



(6.7) 



(6.8) 



(6.9) 



According to the above two relations, the equation of motion (6.3) can be re- 
written as 

su-d(Tu) + ^2n A u-d(iJ, A u) = 0. (6.10) 



In this equation, appears the 1-form 



tv a := [i A u 



(6.11) 



which is called the momentum 1-form of particles of species A. It is called momen- 
tum because in the Hamiltonian formulations mentioned in Sec. 5.1, this 1-form 
is the conjugate of the number density current uau. 

Actually, it is the exterior derivative of tv a which appears in Eq. (6.10): 



d7T 



(6.12) 



This 2-form is called the vorticity 2-form of particles of species A. With this 
definition, Eq. (6.10) becomes 



A 



n A u ■ w A + su ■ d(Tu) = 



(6.13) 



This is the Carter- Lichnerowicz form of the equation of motion for a multi- 
constituent perfect fluid. It has been considered by Lichnerowicz (1967) in the 
case of a single-constituent fluid (N ~ 1) and generalized by Carter (1979, 1989) 
to the multi-constituent case. Let us stress that this is an equation between 1- 
forms. For instance u w A is the 1-form w A (u, .), i.e. at each point p G A4, this is 
the linear application T p — > R, v i— > w A {it, v) . Since w A is antisymmetric (being a 
2-form), w A (u,u) = 0. Hence the Carter-Lichnerowicz equation (6.13) is clearly 
a non-trivial equation only in the three dimensions orthogonal to u. 
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6.2 Canonical form for a simple fluid 

Let us define a simple fluid as a fluid for which the EOS (4.10) takes the form 



e = e{s,n) , (6.14) 



where n is the baryon number density in the fluid rest frame. The simple fluid 
model is valid in two extreme cases: 

• when the reaction rates between the various particle species are very low: 
the composition of matter is then frozen: all the particle number densities 
can be deduced from the baryon number: ua = Yau, with a fixed species 
fraction Ya; 

• when reaction rates between the various particle species are very high, en- 
suring a complete chemical (nuclear) equilibrium. In this case, all the ua 
are uniquely determined by n and s, via ua = Y^ q (s, n) n. 

A special case of a simple fluid is that of barotropic fluid, for which 

e = e(n). (6.15) 

This subcase is particularly relevant for cold dense matter, as in white dwarfs and 
neutron stars. 

Thanks to Eq. (6.14), a simple fluid behaves as if it contains a single particle 
species: the baryons. All the equations derived previously apply, setting N = 1 
(one species) and A = l. 

Since n is the baryon number density, it must obey the fundamental law of 
baryon conservation: 



V • (nu) = . (6.16) 



That this equation does express the conservation of baryon number should be 
obvious after the discussion in Sec. 5.3, from which it follows that cV • (nu) is 
the number of baryons created per unit volume and unit time in a comoving fluid 
element. 

The projection of V • T = along w, Eq. (5.11) then implies 



V • (su) = . (6.17) 



This means that the evolution of a (isolated) simple fluid is necessarily adiabatic. 
On the other side, the Carter-Lichnerowicz equation (6.13) reduces to 

nu-w h + su-d(Tu) = 0, (6.18) 

where we have used the label b (for baryon) instead of the running letter A: 
w h = d(uu) [Eqs. (6.11) and (6.12)], u being the chemical potentials of baryons: 

,:-(£). «U.) 
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Let us rewrite Eq. (6.18) as 

u- [d(uu) + sd(Tu)] = 0, 
where we have introduced the entropy per baryon: 



s 
n 



(6.20) 



(6.21) 



In view of Eq. (6.20), let us define the fluid momentum per baryon 1-form by 

7T := (fj, + Ts)u (6.22) 
and the fluid vorticity 2- form as its exterior derivative: 

(6.23) 



w := d7r 



Notice that w is not equal to the baryon vorticity: w = w h + d(sTu). Since in 
the present case the thermodynamic identity (5.9) reduces to p = Ts + \in — e, we 
have 



(6.24) 



where h is the enthalpy per baryon. Accordingly the fluid momentum per baryon 
1-form is simply 

1 (6.25) 



tv = hu 



By means of formula (2.27), we can expand the exterior derivative d(sTu) as 

d(sTu) = dsA (Tu) + sd(Tu) = Tds Au + sd(Tu), (6.26) 

where the symbol A stands for the exterior product: for any pair (a, b) of 1-forms, 
a A b is the 2-form defined a,s aAb — a(g)b— b(g> a. Thanks to (6.26), the fluid 
vorticity 2-form, given by Eqs. (6.23) and (6.22), can be written as 

w = d(fiu) +sd{Tu) +Tds Am. (6.27) 

Therefore, we may rewrite (6.20) by letting appear w, to get successively 

u-(w- Tds A u) = 0, 

u w - Tu ■ (ds A u) = 0, 

u ■ w — T [(ds, u)u— (u, u)ds] = 0, 

u-w-T[{V a s)u + ds] = 0. (6.28) 

Now from the baryon number and entropy conservation equations (6.16) and 
(6.17), we get 

" " (6.29) 



V a s = 
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i.e. the entropy per baryon is conserved along the fluid lines. Reporting this 
property in Eq. (6.28) leads to the equation of motion 



u ■ w = Tds 



(6.30) 



This equation was first obtained by Lichnerowicz (1967). In the equivalent form 
(assuming T ^ 0), 



u' ■ w = ds, 



with u' := t=u, 

T ■ 



(6.31) 



it has been called a canonical equation of motion by Carter (1979), who has shown 
that it can be derived from a variational principle. 

Owing to its importance, let us make (6.30) explicit in terms of components 
[cf. (2.26) and (2.27)]: 



d 



(hu a ) 



d 
~dx° 



(hUfj) 



= T 



ds 



(6.32) 



6.3 Isen tropic case (barotropic fluid) 

For an isentropic fluid, s — const. The EOS is then barotropic, i.e. it can be 
cast in the form (6.15). For this reason, the isentropic simple fluid is also called 
a single- constituent fluid. In this case, the gradient ds vanishes and the Carter- 
Lichnerowicz equation of motion (6.30) reduces to 



u ■ w 







(6.33) 



This equation has been first exhibited by Synge (1937). Its simplicity is remark- 
able, especially if we compare it to the equivalent Euler form (5.15). Indeed it 
should be noticed that the assumption of a single-constituent fluid leaves the rel- 
ativistic Euler equation as it is written in (5.15), whereas it leads to the simple 
form (6.33) for the Carter-Lichnerowicz equation of motion. 

In the isentropic case, there is a useful relation between the gradient of pressure 
and that of the enthalpy per baryon. Indeed, from Eq. (6.24), we have de + dp = 
d(nh) — ndh + hdn. Substituting Eq. (4.13) for de yields Tds + fidn + dp = 
ndh + hdn. But T ds = T d(ns) = Tsdn since ds = 0. Using Eq. (6.24) again 
then leads to 

dp = ndh, (6.34) 

or cquivalcntly, 



dp 

s + p 



= dkih. 



(6.35) 



If we come back to the relativistic Euler equation (5.15), the above relation shows 
that in the isentropic case, it can be written as the fluid 4-acceleration being the 
orthogonal projection (with respect to u) of a pure gradient (that of — In ft,): 



a = — d In h — (d In h, u) u. 



(6.36) 
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6.4 Newtonian limit: Crocco equation 

Let us go back to the non isentropic case and consider the Newtonian limit of 
the Carter-Lichnerowicz equation (6.30). For this purpose let us assume that the 
gravitational field is weak and static. It is then always possible to find a coordinate 
system (x a ) = (x° = ct, x l ) such that the metric components take the form 

g af3 dx a dx p = -( 1 + 2 Jf) c2dt2 2 j^) -M xW > ( 6 - 37 ) 

where $ is the Newtonian gravitational potential (solution of A$ = AirGp) and 
fij is the flat metric in the usual 3-dimcnsional Euclidean space. For a weak 
gravitational field (Newtonian limit), |$|/c 2 <C 1. The components of the fluid 
4-velocity are deduced from Eq. (3.2): u a = c~ 1 dx a /dr, r being the fluid proper 
time. Thus (recall that x° = ct) 

u°,u°^j, with u° = ^ and v* := (6.38) 

At the Newtonian limit, the v l, s are of course the components of the fluid velocity 
v with respect to the inertial frame defined by the coordinates (x a ). That the coor- 
dinates (x a ) are inertial in the Newtonian limit is obvious from the form (6.37) of 
the metric, which is clearly Minkowskian when <i> — > 0. Consistent with the New- 
tonian limit, we assume that \v\/c <C 1. The normalization relation g a pu a uP = —1 
along with (6.38) enables us to express u° in terms of $ and v. To the first order 
in <3>/c 2 and v ■ v/c 2 = VjV^ /c 2 ( 4 ), we get 

To that order of approximation, we may set u° ~ 1 in the spatial part of u a and 
rewrite (6.38) as 

The components of u are obtained from u a — gapu 13 , with g a p given by (6.37). 
One gets 



U a ~ (u , ^ ~ ^~ 



To form the fluid vorticity w we need the enthalpy per baryon h. By combining 
Eq. (6.24) with Eq. (4.11) written as e = rrihnc 2 +e[ n t (where Wb is the mean mass 
of one baryon: mb — 1.66 x 10 -27 kg), we get 

h = m, h c 2 ( 1 + \ j , (6.42) 



4 the indices of v' are lowered by the flat metric: Vi := fijvi 
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where H is the non-relativistic (i.e. excluding the rest-mass energy) specific en- 
thalpy (i.e. enthalpy per unit mass): 



H := 



From (6.32), we have, for i 6 {1,2,3}, 



£int + V 



(6.43) 



fit 



= u 



l c§-t ihUl) "^ ikUa) 



Plugging Eqs. (6.40), (6.41) and (6.42) yields 



Hi 



d_ 

Wt 

\- 



1 + -7T \Vi 



(6.44) 



1 H \ 



d 

dx 1 



l + -n)Vi 



(6.45) 



At the Newtonian limit, the terms u° and H/c 2 in the above equation can be set 
to respectively 1 and 0. Moreover, thanks to (6.41), 



_d_ 

dx 



- [(c 2 +H)u Q ] 



$ VjV 

(c ; +Jff)( l + _ + ^_ 



d 

dx 1 



Finally we get 



d 



dt dx 1 



——— -+- — — r [ H -+- -Viv> + $ + v 3 -4 - -4 



da;-? 9a;' 



(6.46) 



(6.47) 



The last term can be expressed in terms of the cross product between v and its 
the (3-dimensional) curl: 



dvi dvj 
dxi dx 1 



(v x cur\v) i . 



(6.48) 



In view of (6.47) and (6.48), we conclude that the Newtonian limit of the Carter- 
Lichncrowicz canonical equation (6.30) is 



dvi d 



~dt + dx~^ H+ - VlV ° + * ' ~ {tJ X CnTl ^ = 



d~s 



dx 1 



(6.49) 



where s := s/m b is the specific entropy (i.e. entropy per unit mass). Equation 
(6.49) is known as the Crocco equation [see e.g. (Rieutord 1997)]. It is of course 
an alternative form of the classical Euler equation in the gravitational potential 
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7 Conservation theorems 

In this section, we illustrate the power of the Carter-Lichnerowicz equation by 
deriving from it various conservation laws in a very easy way. We consider a 
simple fluid, i.e. the EOS depends only on the baryon number density and the 
entropy density [Eq. (6.14)]. 

7.1 Relativistic Bernoulli theorem 

7.1.1 Conserved quantity associated with a spacetime symmetry 

Let us suppose that the spacetime (M,g) has some symmetry described by the 
invariance under the action of a one-parameter group Q: for instance Q = (R, +) 
for stationarity (invariance by translation along timelike curves) or Q = SO (2) for 
axisymmetry (invariance by rotation around some axis). Then one can associate 
to Q a vector field £ such that an infinitesimal transformation of parameter e in the 
group Q corresponds to the infinitesimal displacement e£. In particular the field 
lines of £ are the trajectories (also called orbits) of Q. £ is called a generator of 
the symmetry group Q or a Killing vector of spacetime. That the metric tensor g 
remains invariant under Q is then expressed by the vanishing of the Lie derivative 
of g along £: 

(7-1) 



C e g = 



Expressing the Lie derivative via Eq. (2.15) with the partial derivatives replaced 
by covariant ones [cf. remark below Eq. (6.4)], we immediately get that (7.1) is 
equivalent to the following requirement on the 1-form £ associated to | by the 
metric duality: 

V a & + = 0. (7.2) 

Equation (7.2) is called the Killing equation. It fully characterizes Killing vectors 
in a given spacetime. 

The invariance of the fluid under the symmetry group Q amounts to the van- 
ishing of the Lie derivative along | of all the tensor fields associated with matter. 
In particular, for the fluid momentum per baryon 1-form 7r introduced in Sec. 6.2: 

Cgir = 0. (7.3) 

By means of the Cartan identity (2.32), this equation is recast as 

fw + d(n,i) = 0, (7.4) 

where we have replaced the exterior derivative dir by the vorticity 2-form w [cf. 
Eq. (6.23)] and we have written £ -it = (ir, (scalar field resulting from the action 
of the 1-form 7r on the vector £). The left-hand side of Eq. (7.4) is a 1-form. Let 
us apply it to the vector u: 

to(£tI) + Vfi(7r,3=0 (7.5) 
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Now, since w is antisymmetric, w(£,u) = —w(u,£) and we may use the Carter- 
Lichnerowicz equation of motion (6.30) which involves w(u, .) = u ■ w to get 

-T(ds,i) + V a (TT,i)=0. (7.6) 

But (ds, £) = Cgs and, by the fluid symmetry under Q, Cgs = 0. Therefore there 
remains 

Va(7r,€) = 0, (7.7) 
which, thanks to Eq. (6.25), we may rewrite as 



V fi (/i£-u) = 



(7.8) 



We thus have established that if £ is a symmetry generator of spacetime, the scalar 
field h£ • u remains constant along the flow lines. 

The reader with a basic knowledge of relativity must have noticed the simi- 
larity with the existence of conserved quantities along the geodesies in symmetric 
spacetimes: if £ is a Killing vector, it is well known that the quantity £ • u is 
conserved along any timelike geodesic (u being the 4-velocity associated with the 
geodesic) [see e.g. Chap. 8 of (Hartle 2003)]. In the present case, it is not the 
quantity £ • u which is conserved along the flow lines but h£ ■ u. The "correcting 
factor" h arises because the fluid worldlines are not geodesies due to the pressure 
in the fluid. As shown by the relativistic Euler equation (5.15), they are geodesies 
(a = 0) only if p is constant (for instance p = 0) . 



7.1.2 Stationary case: relativistic Bernoulli theorem 

In the case where the Killing vector £ is timelike, the spacetime is said to be 
stationary and Eq. (7.8) constitutes the relativistic generalization of the classi- 
cal Bernoulli theorem. It was first established by Lichnerowicz (1940) (see also 
Lichnerowicz 1941), the special relativistic subcase (flat spacetime) having been 
obtained previously by Synge (1937). 

By means of the formulae established in Sec. 6.4, it is easy to see that at 
the Newtonian limit, Eq. (7.8) does reduce to the well-known Bernoulli theorem. 
Indeed, considering the coordinate system (x a ) given by Eq. (6.37), the Killing 
vector £ corresponds to the invariance by translation in the t direction, so that we 
have £ = d/dx° = c~ 1 d/dt. The components of £ with respect to the coordinates 
(x a ) are thus simply 

r = (1,0,0,0). (7.9) 



Accordingly 



hi- u - hu a C = hu ~ -m h c 2 + - l + - + (7.10) 



H\ ( _ $ VjV J 
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where we have used Eq. (6.42) for h and Eq. (6.41) for uq. Expanding (7.10) to 
first order in c~ 2 , we get 

h£-u~ -m b (c 2 + H + ^-VjV j + $ ) . (7.11) 



2 



Since thanks to Eq. (6.40), 



v a (/»€ • «) = • «) = - ^(^ • «) • «)- ( 7 - 12 ) 

S v ' 

=0 

we conclude that the Newtonian limit of Eq (7.8) is 

dx l V 2 

i.e. we recover the classical Bernoulli theorem for a stationary flow 



v*— ( H + -v 3 v j + $ ) = 0, (7.13) 



7.1.3 Axisymmetric flow 

In the case where the spacetime is axisymmetric (but not necessarily stationary), 
there exists a coordinate system of spherical type x a — (x° = ct, r, 6, ip) such that 
the Killing vector is 

* = h (7 - 14) 

The conserved quantity h £ • u is then interpretable as the angular momentum per 
baryon. Indeed its Newtonian limit is 

u = hu a £ a = hu v ~ m b c 2 ^1 + ^ ~ m h cv v , (7.15) 

where we have used Eq. (6.41) to replace u v by v v /c. In terms of the components 
of the fluid velocity in an orthonormal frame, one has v v = rsinOv^, so that 

■ u — c x r smOm^v^y (7.16) 

Hence, up to a factor c, the conserved quantity is the z-component of the angular 
momentum of one baryon. 

7.2 Irrotational flow 

A simple fluid is said to be irrotational iff its vorticity 2-form vanishes identically: 

(7.17) 



w = 



It is easy to see that this implies the vanishing of the kinematical vorticity vector 
lj defined by 

^ := \^u,V p u a = \^u^ p . (7.18) 
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In this formula, e apf " T stands for the components of the alternating type-(4, 0) 
tensor e that is related to the volume element 4-form e associated with g by 
e af3 ^ 5 e a p~ ( 8 = —4!. Equivalently, e is such that for any basis of 1-forms (e Q ) dual 
to a right-handed orthonormal vector basis (e a ), then e(e°, e 1 , e 2 , e 3 ) = 1. Notice 
that the second equality in Eq. (7.18) results from the antisymmetry of e combined 
with the symmetry of the Christoffcl symbols in their lower indices. From the 
alternating character of e, the kinematical vorticity vector u5 is by construction 
orthogonal to the 4-velocity: 

u-w = 0. (7.19) 

Moreover, at the non-relativistic limit, cD is nothing but the curl of the fluid ve- 
locity: 

u~-cui\v. (7.20) 



c 

That w = implies u = 0, as stated above, results from the relation 

1 



Ah 



e a ^ p(7 u u u 



jj, ^ pa •) 



(7.21) 



which is an easy consequence of w pa = d{hu a )/ dx p — d{hu a )/ dx a [Eq. (6.23)]. 
From a geometrical point of view, the vanishing of c3 implies that the fluid world- 
lines are orthogonal to a family of (spacelike) hypersurfaces (submanifolds of M. 
of dimension 3). 

The vanishing of the vorticity 2-form for an irrotational fluid, Eq. (7.17), means 
that that the fluid momentum per baryon 1-form 7T is closed: dir = 0. By Poincare 
lemma (cf. Sec. 2.5), there exists then a scalar field \P such that 



7T = 



i.e. 



hu = d* 



(7.22) 



The scalar field iff is called the potential of the flow. Notice the difference with 
the Newtonian case: a relativistic irrotational flow is such that h u is a gradient, 
not u alone. Of course at the Newtonian limit h — > m^c 2 = const, so that the two 
properties coincide. 

For an irrotational fluid, the Carter-Lichnerowicz equation of motion (6.30) 
reduces to 

Tds = 0. (7.23) 

Hence the fluid must either have a zero temperature or be isentropic. The con- 
straint on ^ arises from the baryon number conservation, Eq. (6.16). Indeed, we 
deduce from Eq. (7.22) that 

M = iv*, (7.24) 

where VvP denotes the vector associated to the gradient 1-form d^ = V* by the 
standard metric duality, the components of being V a $ = g ap V^. Inserting 
(7.24) into the baryon number conservation equation (6.16) yields 

^□* + VQ) • V* = 0, (7.25) 
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Velocity w.r.t corotating frame (z = 0) 




Fig. 5. Velocity with respect to a co-orbiting observer for irrotational binary relativistic 
stars. The figure is drawn in the orbital plane. The velocity field has been obtained 
by numerically solving Eq. (7.25) for the fluid potential ty. The first integral of motion 
(7.27) provided by the helical Killing vector has been used to get the enthalpy per baryon 
h. The density profile in the stars is then deduced from the EOS [from (Gourgoulhon et 
al. 2001)]. 



where □ is the d'Alembertian operator associated with the metric g: □ := V- V = 

Let us now suppose that the spacetime possesses some symmetry described by 
the Killing vector £. Then Eq. (7.4) applies. Since w = in the present case, it 
reduces to 

d<7r,O = 0. (7.26) 
We conclude that the scalar field (7r,£) is constant, or equivalently 

h£-u = const. (7.27) 



Hence for an irrotational flow, the quantity h£ ■ u is a global constant, and not 
merely a constant along each fluid line which may vary from a fluid line to another 
one. One says that h £ ■ u is a first integral of motion. This property of irrotational 
relativistic fluids was first established by Lichnerowicz (1941), the special relativis- 
tic subcase (flat spacetime) having been proved previously by Synge (1937). 

As an illustration of the use of the integral of motion (7.27), let us consider 
the problem of equilibrium configurations of irrotational relativistic stars in binary 
systems. This problem is particularly relevant for describing the last stages of the 
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slow inspiral of binary neutron stars, which are expected to be one of the strongest 
sources of gravitational waves for the interferometric detectors LIGO, GEO600, 
and VIRGO (sec Baumgartc & Shapiro (2003) for a review about relativistic binary 
systems). Indeed the shear viscosity of nuclear matter is not sufficient to synchro- 
nize the spin of each star with the orbital motion within the short timescale of 
the gravitational radiation-driven inspiral. Therefore contrary to ordinary stars, 
close binary system of neutron stars are not in synchronized rotation. Rather if 
the spin frequency of each neutron star is initially low (typically 1 Hz), the orbital 
frequency in the last stages is so high (in the kHz regime), that it is a good approxi- 
mation to consider that the fluid in each star is irrotational. Besides, the spacetimc 
containing an orbiting binary system has a priori no symmetry, due to the emis- 
sion of gravitational wave. However, in the inspiral phase, one may approximate 
the evolution of the system by a sequence of equilibrium configurations consisting 
of exactly circular orbits. Such a configuration possesses a Killing vector, which 
is helical, being of the type £ = d/dt + Sld/dip, where Q is the orbital angular 
velocity (see Friedman et al. 2002 for more details). This Killing vector provides 
the first integral of motion (7.27) which permits to solve the problem (cf. Fig. 5). 
It is worth noticing that the derivation of the first integral of motion directly from 
the relativistic Euler equation (5.15), i.e. without using the Carter-Lichnerowicz 
equation, is quite lengthy (Shibata 1998, Teukolsky 1998). 



7.3 Rigid motion 

Another interesting case in which there exists a first integral of motion is that of 
a rigid isentropic flow. We say that the fluid is in rigid motion iff (i) there exists 
a Killing vector £ and (ii) the fluid 4-velocity is collinear to that vector: 



u = A£, (7.28) 



where A is some scalar field (not assumed to be constant). Notice that this relation 
implies that the Killing vector £ is timelike in the region occupied by the fluid. 
Moreover, the normalization relation u ■ u = —1 implies that A is related to the 
scalar square of £ via 

/ -, -A 1/2 

A = (-£ • • (7.29) 

The denomination rigid stems from the fact that (7.28) in conjunction with the 
Killing equation (7.2) implies 5 that both the expansion rate := V • u [cf. 
Eq. (5.13)] and the shear tensor (see e.g. Ehlers 1961) 



1,_ _ ™, 1 

2 



(Tap ■= o ( + Vv^P^P^ - -eP af3 (7.30) 



vanish identically for such a fluid. 



5 the reverse has been proved to hold for an isentropic fluid (Salzman & Taub 1954) 
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Equation (7.4) along with Eq. (7.28) results in 



\u-w + d(n,O=0. (7.31) 
A 

Then, the Carter-Lichnerowicz equation (6.30) yields 

^ds + d(7r,£) = 0. (7.32) 
A 

If we assume that the fluid is isentropic, ds = and we get the same first integral 
of motion than in the irrotational case: 



h£ ■ u — const. 



(7.33) 



This first integral of motion has been massively used to compute stationary and 
axisymmetric configurations of rotating stars in general relativity (see Stergioulas 
2003 for a review). In this case, the Killing vector £ is 

£ = ^station ~\~ ^£axisym> (7.34) 

where ^station and £ a xisym are the Killing vectors associated with respectively sta- 
tionarity and axisymmetry. Note that the isentropic assumption is excellent for 
neutron stars which are cold objects. 

7.4 First integral of motion in symmetric spacetimes 

The irrotational motion in presence of a Killing vector £ and the isentropic rigid 
motion treated above are actually subcases of flows that satisfy the condition 



£ • w = , (7.35) 



which is necessary and sufficient for (n, £) = h £ • u to be a first integral of motion. 
This property follows immediately from Eq. (7.4) (i.e. the symmetry property 
Cg-n = re-expressed via the Cartan identity) and was first noticed by Lich- 
nerowicz (1955). For an irrotational motion, Eq. (7.35) holds trivially because 
w = 0, whereas for an isentropic rigid motion it holds thanks to the isentropic 
Carter-Lichnerowicz equation of motion (6.33) with u = A£. 



7.5 Relativistic Kelvin theorem 

Here we do no longer suppose that the spacetime has any symmetry. The only 
restriction that we set is that the fluid must be isentropic, as discussed in Sec. 6.3. 
The Carter-Lichnerowicz equation of motion (6.33) leads then very easily to a 
relativistic generalization of Kelvin theorem about conservation of circulation. In- 
deed, if we apply Cartan identity (2.32) to express the Lie derivative of the fluid 
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vorticity 2-form w along the vector au (where a is any non- vanishing scalar field) , 
we get 

C a n w = au • dw +d(auji)). (7.36) 

The first "= 0" results from dw = dd7r = [nilpotent character of the exterior 
derivative, cf. Eq. (2.30)], whereas the second "= 0" is the isentropic Carter- 
Lichnerowicz equation (6.33). Hence 



C na w = 



(7.37) 



This constitutes a relativistic generalization of Hclmholtz's vorticity equation 

^ =C5- Vv-v- Vw-(divtT)w, (7.38) 

which governs the evolution of uj :— curliT [cf. Eq. (7.20)]. 

The fluid circulation around a closed curve C C M is defined as the integral 
of the fluid momentum per baryon along C: 



C(C) := / tt 
Jc 



(7.39) 



Let us recall that C being a 1-dimensional manifold and tv a 1-form the above 
integral is well defined, independently of any length element on C. However to 
make the link with traditional notations in classical hydrodynamics, we may write 
tv = hu [Eq. (6.25)] and let appear the vector u (4- velocity) associated to the 
1-form u by the standard metric duality. Hence we can rewrite (7.39) as 

C(C) = [ hu-dl. (7.40) 
Jc 

This writing makes an explicit use of the metric tensor g (in the scalar product 
between u and the small displacement d£). 

Let S be a (2-dimensional) compact surface the boundary of which is C: C = 
dS. Then by the Stokes theorem (2.31) and the definition of w, 

C(C) = [ dTr = / w. (7.41) 
Js Js 

We consider now that the loop C is dragged along the fluid worldlines. This 
means that we consider a 1-parameter family of loops C(A) that is generated from 
a initial loop C(0) nowhere tangent to u by displacing each point of C(0) by some 
distance along the field lines of u. We consider as well a family of surfaces S(X) 
such that dS(X) = C(A). We can parametrize each fluid worldline that is cut by 
5(A) by the parameter A instead of the proper time r. The corresponding tangent 
vector is then v = au, where a := dr/dX (the derivative being taken along a 
given fluid worldline). From the very definition of the Lie derivative (cf. Sec. 2.4.1 
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where the Lie derivative of a vector has been defined from the dragging of the 
vector along the flow lines), 




(7.42) 



From Eq. (7.37), we conclude 



(7.43) 



This is the relativistic Kelvin theorem. It is very easy to show that in the Newtonian 
limit it reduces to the classical Kelvin theorem. Indeed, choosing A = r, the non- 
relativistic limit yields t = t, where t is the absolute time of Newtonian physics. 
Then each curve C(t) lies in the hypersurface t = const (the "space" at the instant 
t, cf. Sec. 2.1). Consequently, the scalar product u ■ dl in (7.40) involves only the 
spatial components of u, which according to Eq. (6.40) are u l ~ v l /c. Moreover 
the Newtonian limit of h is nibc 2 [cf. Eq. (6.42)], so that (7.40) becomes 

C(C) ~ m h c ( v-dt (7.44) 
Jc 

Up to the constant factor m^c we recognize the classical expression for the fluid 
circulation around the circuit C. Equation (7.43) reduces than to the classical 
Kelvin theorem expressing the constancy of the fluid circulation around a closed 
loop which is comoving with the fluid. 

7.6 Other conservation laws 

The Carter-Lichnerowicz equation enables one to get easily other relativistic con- 
servation laws, such as the conservation of helicity or the conservation of enstrophy. 
We shall not discuss them in this introductory lecture and refer the reader to ar- 
ticles by Carter (1979, 1989), Katz (1984) or Bckenstcin (1987). 

8 Conclusions 

The Carter-Lichnerowicz formulation is well adapted to a first course in relativistic 
hydrodynamics. Among other things, it uses a clear separation between what is 
a vector and what is a 1-form, which has a deep physical significance (as could 
also be seen from the variational formulations of hydrodynamics mentioned in 
Sec. 5.1). For instance, velocities are fundamentally vectors, whereas momenta 
are fundamentally 1-forms. On the contrary, the "standard" tensor calculus mixes 
very often the concepts of vector and 1-form, via an immoderate use of the metric 
tensor. Moreover, we hope that the reader is now convinced that the Carter- 
Lichnerowicz approach greatly facilitates the derivation of conservation laws. It 
must also be said that, although we have not discussed it here, this formulation 
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can be applied directly to non-relativistic hydrodynamics, by introducing exterior 
calculus on the Newtonian spacetime, and turns out to be very fruitful (Carter 
& Gaffet 1988; Prix 2004; Carter & Chamel 2004; Chamel 2004). Besides, it is 
worth to mention that the Carter-Lichnerowicz approach can also be extended to 
relativistic magnetohydrodynamics (Lichnerowicz 1967 and Sec. 9 of Carter et al. 
2006). 

In this introductory lecture, we have omitted important topics, among which 
relativistic shock waves (see e.g. Marti & Miiller 2003; Font 2003, Anile 1989), 
instabilities in rotating relativistic fluids (see e.g. Stergioulas 2003; Andersson 
2003; Villain 2006), and superfluidity (see e.g. Carter & Langlois 1998, Prix et 
al. 2005). Also we have not discussed much astrophysical applications. We may 
refer the interested reader to the review article by Font (2003) for relativistic hy- 
drodynamics in strong gravitational fields, to Shibata et al. (2005) for some recent 
application to the merger of binary neutron stars, and to Baiotti et al. (2005), 
Dimmelmeier et al. (2005), and Shibata & Sekiguchi (2005) for applications to 
gravitational collapse. Regarding the treatment of relativistic jets, which requires 
only special relativity, we may mention Sauty et al. (2004) and Marti & Miiller 
(2003) for reviews of respectively analytical and numerical approaches, as well as 
Alloy & Rezzola (2006) for an example of recent work. 

It is a pleasure to thank Berangcre Dubrulle and Michel Rieutord for having organized the very- 
successful Cargese school on Astrophysical Fluid Dynamics. I warmly thank Brandon Carter 
for fruitful discussions and for reading the manuscript. I am extremely grateful to Silvano 
Bonazzola for having introduced me to relativistic hydrodynamics (among other topics !), and 
for his constant stimulation and inspiration. In the spirit of the Cargese school, I dedicate this 
article to him, recalling that his very first scientific paper (Bonazzola 1962) regarded the link 
between physical measurements and geometrical operations in spacetime, like the orthogonal 
decomposition of the 4-velocity which we discussed in Sec. 3. 
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